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Abstract In this paper we report a computational study of the effects of 
strain on the conductivity of Y-doped ceria (YDC). This material was chosen 
as it is of technological interest in the field of Solid Oxide Fuel Cells (SOFCs). 
— i The simulations were performed under realistic operational temperatures and 

strain (e) levels. For bulk and thin film YDC, the results show that tensile 
strain leads to conductivity enhancements of up to 3.5 x and 1.44 x, respec- 
tively. The magnitude of these enhancements is in agreement with recent ex- 
perimental and computational evidence. In addition, the methods presented 



< 



M. Burbano 

School of Chemistry and CRANN 
. t Trinity College Dublin 

Dublin 2 

S Ireland 

Tel.: +353-877698113 
——t E-mail: burbanom@tcd.ie 

£^ D. Marrocchelli 

School of Chemistry and CRANN 
(^J Trinity College Dublin 

I— I Dublin 2 

Ireland 

Department of Materials Science and Engineering 
Department of Nuclear Science and Engineering 
Massachusetts Institute of Technology 

f<~) Cambridge, MA 

^- USA 

VP E-mail: dmarrocc@mit.edu 

^-" G.W. Watson 

School of Chemistry and CRANN 
/■^-\ Trinity College Dublin 

Dublin 2 
• • Ireland 

> Tel.: +353-18961357 

E-mail: watsong@tcd.ie 

X 



Mario Burbano ct al. 



herein allowed us to identify enhanced ionic conductivity in the surface regions 
of YDC slabs and its anisotropic character. 



Keywords SOFC • ceria • strain • MD 



1 Introduction 



Solid Oxide Fuel Cells (SOFCs) are electrochemical devices which efficiently 
generate electricity from a variety of fuels pQ . As such, they have the potential 
to play a key role in the energy production landscape in years to come 0[3J. 
Significant research efforts have been focused, over the past decades, on im- 
proving the performance of the oxide ion conductors used either as electrolytes 
or cathodes in these devices [TJIMS]. Operating SOFCs at temperatures lower 
than those required by current Yttria Stabilized Zirconia (YSZ) technologies 
(~ 1273 K) would allow for the use of cheaper materials in the fabrication of 
the cells, as well as, facilitating their possible use in portable electronics, e.g. \i- 
SOFCs 6J. In particular, Rare Earth doped ceria (REDC), Cei_ a: RE 2: 02- 2 ;/2, 
where RE = Y (YDC), La (LDC), Nd (NDC), Sm (SDC), Gd (GDC), etc has 
been studied extensively given its relatively high ionic conductivity (ctj) in the 
intermediate temperature range (873 K - 1073 K). This process is known to 
occur by means of vacancy migration in the anion sublattice 7,8. Equation 
[l] illustrates, in Kroger- Vink notation, the formation of charge compensating 
vacancies (Vq) upon substitution of a host Ce + cation by a dopant RE 3+ 
cation (REq ). 

RE 2 3 + (Ce0 2 ) n -> (Ce0 2 ) n _ 2 + 2REc c + 30£ + V (1) 

The interactions between these defects can be classified as cation-vacancy, 
vacancy- vacancy and cation-cation p iTUIITTllT2llT5llTiIITMT()] . In practice, these 
interactions can prove deleterious to the ionic conductivity because they de- 
crease the number of mobile vacancies available, hence the need for high op- 
erating temperatures |17) . 

From its outset, much of the research within the context of electrolyte ap- 
plications has focused on the optimization of ionic conductivity from a com- 
positional point of view, i.e. the identification of the dopant species and the 
concentration at which defect interactions are minimized, especially those be- 
tween dopants and vacancies [10, 18, 19, 20, 21,22.. In the case of ceria, these 
studies have identified RE elements, such as Gd 3+ , Y 3+ , Sm 3+ and Pm 3+ as 
the best candidate dopants, given that their radius mismatch with the host 
cation (Ce 4+ ) balances the competing electrostatic and elastic components of 
the defect interactions which control their association [53]. Nevertheless, re- 
cent studies have shown that in the limit where cation- vacancy interactions are 
reduced to a minimum, it is vacancy-vacancy association which ultimately de- 
termines the ionic conductivity drop as a function of dopant concentration in 
fluorite-structured materials, such as, YDC [16], YSZ and Scandia Stabilized 
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Zirconia (ScSZ) |15) . This means that further improvements in the perfor- 
mance of IT-SOFC electrolytes must come from the exploitation of different 
optimization routes [2KM3]. 

Current thin film deposition techniques allow for the fabrication of epi- 
taxially grown oxides with high levels of control over the microstructure, sto- 
ichiometry and lattice mismatch of individual layers 24,25,26,27,28 . These 
developments have attracted much attention to this area of research because 
they might open a new avenue for SOFC electrolyte optimization in a way sim- 
ilar to that which has been achieved in Si-based semiconductors [29] . Changes 
in the ionic conductivity of thin film electrolytes, whether detrimental or ben- 
eficial, are ascribed to interfacial effects present at phase boundaries. One 
such effect is the formation of space charge regions which have been shown 
to cause a dramatic increase in the fluoride ion conductivity of BaF2/CaF2 
heterostructures [30,31 . However, the the Debye screening length of extrin- 
sic ionic conductors, such as, doped ceria/zirconia, is too short (given that it 
is inversely proportional to the square root of the defect concentration) and, 
as a result, space charge regions are not expected to play an important role 
in these systems 24,32,28 . Coherent growth of a thin film on a substrate, 
or heterostructure, leads to strain (e) due to the structural mismatch which 
arises from differences in their lattice vectors (Equation [2]). 

6 =*^° (2) 

a 

Thin films can elastically accommodate strains of e ~ 0.03, depending on the 
thickness of the film and the elastic properties of the material [2"5"1[3"3"] . Strain 
levels beyond this value are typically released by the films through the forma- 
tion of dislocation networks, defect clustering at the interface and rotations 
with respect to the substrate or other layers [3~4"l[3"5"ll3"6"] . 

Substantial research efforts have been devoted, over the last decade, to the 
elucidation of how interfacial effects alter ionic conductivity. Many such stud- 
ies have focused on YSZ as it is the prevalent electrolyte in SOFCs. The results 
reported in the literature for this material show significant scatter in the range 
of conductivities achieved through the formation of interfaces. Experimental 
investigations have observed conductivity enhancements varying from a colos- 
sal 10 8 increase in ctj [37] down to neglegible changes [38] and myriad values 
in between P3 l l59 l l52 l [2^ 1 l4l? l l4T ] H2 l li3 l li4] . Early theoretical work using Molec- 
ular Dynamics (MD) simulations found only a small enhancement in Oi for 
YSZ mediated by a lowering of the activation barrier |45j , while more recent 
Density Functional Theory (DFT) calculations have predicted a maximum en- 
hancement of up to four orders of magnitude in this property, at 400 K for 
relatively high strains (e = 4%) [4"6"ll4"7] . 

In the case of REDC, the conductivity enhancements observed have been 
more modest. For example, Chen et al. studied thin GDC films grown on MgO 
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(lattice mismatch of 28%) and found only a small increase in the ionic con- 
ductivity with respect to the bulk system [34]. Suzuki et al. [48] spin coated 
sapphire substrates with CeC>2 thin films, as well as, GDC films of different 
concentrations. Their study found that the conductivity increased with dimin- 
ishing film thickness. However, this change may be attributed to a lowering 
of the oxygen vacancy formation energy which leads to reduction of Ce 4+ to 
Ce 3+ (Equation pi), and thus, to electronic conductivity [49j . 

2Ce£ c + Og o V + \ 2 + 2Cec c (3) 

This is an undesirable effect in SOFC electrolytes as it causes an internal short 
circuit and cell delamination through lattice expansion [50,25,51,52.53,54, 
155] . Similar results were obtained by Perkins et at, who saw the formation of 
discrete micro-domains containing Ce 3+ in SDC/Ce02 heterostructures (e ~ 
0.0035) grown on MgO substrates [55]. Studies of the ionic conductivity in 
epitaxial SDC films grown on MgO substrates using SrTi03 buffer layers (e ~ 
0.016) found a conductivity of 0.07 S cm" 1 at 973 K compared to ~0.02 S cm" 1 
for dense polycrystalline pellets at the same temperature [27 . In the case of 
STO-buffered SDC/YSZ heterostructures (e - 0.055)) grown on MgO, the 
same group found a conductivity increase of two orders of magnitude with re- 
spect to SDC polycrystaline pellets and about one order of magnitude increase 
compared to either SDC or YSZ thin films [2H]. A recent computational study, 
which employed static calculations with interatomic potentials, predicted an 
enhancement in Oi of up to four orders of magnitude when tensile strain (e = 
0.040) was applied to bulk Ce0 2 [57] , 

This brief survey of our current understanding of interfacial effects high- 
lights the need for further investigation in this area of materials research. To 
this end, computer simulations can play an important role because the vari- 
ous processes which underlie complex phenomena can be treated in a direct 
and controlled manner, thus, making it possible to evaluate their contribu- 
tion to the overall changes observed in experiments. In this paper we have 
assessed how the ionic conductivity in YDC is modified when bulk and thin 
films are subjected to isotropic and anisotropic biaxial strain, respectively. Our 
study has a series of distinctive features that set it apart from previous work. 
Firstly, we use molecular dynamics simulations coupled with accurate dipole- 
polarizable interatomic potentials derived directly from ab initio calculations 
[55UTB] . Secondly, our use of the slab method allows us to account for the re- 
laxation perpendicular to the plane where strain is applied (Poisson effect), as 
well as the effect on the ionic conductivity that results from the presence of 
surfaces on this material; all under realistic dopant concentrations (10-20 %) 
and high temperatures (1273-1673 K). Also, this paper differs from previous 
ones [J61I47U571I55] in that we only study small strain levels (e < 0.021), that 
can be elastically accommodated in thin films. Finally, in our simulations, the 
Ce cation reduction reaction (Ce 4+ to Ce 3+ ), observed in some studies [48], is 
not allowed. This makes it possible to isolate the effects of strain on the ionic 
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conductivity only. 

2 Methods 

2.1 Interatomic Potential 

The computational methods used in this study are well established and have 
been described elsewhere for doped ceria and other oxides [60,58,16,61,62], as 
well as fluoride systems [5511 64 j . The interatomic potential that was employed 
is known as the DIPole Polarizable Ion Model (DIPPIM) [65 . In this model 
the various ionic species are assigned their formal valence charges (Ce 4+ , Y 3+ , 
and 2 ~ ) with the inclusion of the polarization effects that result from the in- 
duction of dipoles on the ions. The DIPPIM parameters were derived directly 
from hybrid density functional theory (h-DFT) calculations using the Heyd, 
Scuzeria, Ernzerhof (HSE06) functional [B6ll67| . as implemented in the VASP 
code 68J. These static h-DFT calculations provided forces and dipoles of each 
ion, and were used to optimize the h-DIPPIM (HSE-DIPPIM) parameters. 
The procedure employed to generate a set of potentials for Ce02-Zr02 solid 
solutions doped (or reduced) with trivalent RE cations (Sc 3+ , Y 3+ , Gd 3+ , 
Sm 3 +, Nd 3 +, Ce 3+ , La 3+ ) will be the subject of a future publication [55 . The 
inclusion of a fraction of nonlocal Hartree-Fock exchange to standard DFT in 
functionals such as HSE is known to be necessary to correctly describe the 
electronic strucuture of lanthanide oxides such as ceria [7D1I7TJ . This is im- 
portant in cases where highly correlated /-electrons are present. However, for 
YDC this potential represents an improvement in terms of a closer agreement 
to experimental lattice constants [70]. Table [I] presents the parameters for the 
h-DIPPIM interatomic potential used in this study. 



2.2 MD Simulation Details 

All Molecular Dynamics simulations were performed with an in-house code 
(PIMAIM) which uses three-dimensional periodic boundary conditions. The 
ionic conductivities of three different dopant concentrations were studied, 
namely Cei_ a; Y :E 02_ a ;/2 where x = 0.08,0.12,0.18. This choice of values origi- 
nated from our previous work, where it was found that the conductivity max- 
ium in YDC is at x = 0.12 [16]. The other two values for the dopant concen- 
tration were included in order to investigate if straining YDC causes a shift in 
the position of the conductivity maximum. As was the case in previous stud- 
ies, the dopant cations and their corresponding charge compensating vacan- 
cies were distributed at random within their respective sublattices. Short MD 
simulations were carried out for each concentration at constant temperature 
and pressure (NPT ensemble) |72j in order to obtain the equilibrium lattice 
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Table 1 Parameters for the h-DIPPIM potential. All values are in atomic units. 



constants at three temperatures, T = 1673 K, 1473 K and 1273 K. Isotropic 
tensile strain was simulated by expanding the axes of 6 x 6 x 6 fluorite super- 
cells (henceforth referred to as bulk) under constant temperature and volume 
conditions (NVT ensemble). Each bulk simulation cell was subjected to three 
different strain levels e x — e y — e z — {0.007, 0.014, 0.021} for each tempera- 
ture. YDC thin films of the same concentrations as bulk were simulated using 
slabs that exposed the (111) surface which has been found to be the most 
stable for ceria both in experiments and simulations [73 , 74 . The size of the 
simulation cells along x and y were ~ 30 A and ~ 26 A (Figure [lp), respectively, 
depending on the dopant concentration and temperature. The length along the 
z-axis for the slab simulation cells was 80 A. This corresponded to a vacuum 
gap of ~ 40 A depending on the dopant concentration and temperature (Pic- 
ture lib) . The vacuum gap is necessary in order to avoid interactions between 
periodic images of the slabs along this direction. By using slabs with these 
dimensions it was possible to have similar numbers of ions of each species 
as in the bulk YDC cells. Table [2] reports the number of atoms from each 
species in these systems. The simulation of anisotropic biaxial strain on YDC 
was performed by applying every combination of strain to the slabs, where 
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ex,e y = {0.007,0.014,0.021} (Picture [TJi) . These simulations were run in the 
NVT ensemble resulting in the lattice vectors along the x and y directions of 
the slabs being fixed while relaxation was allowed along the direction normal 
to the slab surface (/-axis) in order to account for the Poisson effect (Picture 
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Fig. 1 a) CeC>2 slab xj/-plane b) Slab yz-plane 
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x in Cei-zY^Oa-^ 
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1631 


0.12 
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1677 


734 


98 


1615 


0.18 


710 


154 


1651 


684 


148 


1590 



Table 2 Number of ions from each species in the bulk and slab simulation cells for each 
concentration of YDC. 



The use of the DIPPIM interatomic potential made it possible to perform 
the long simulations required to measure the diffusion coefficients of these sys- 
tems, which in the case of the lowest temperature considered was 600 ps with a 
timestep of 1 fs for all calculations. Coulombic and dispersion interactions were 
handled using Ewald summations, while the short-range part of the potential 
was truncated at 12.7 A. 
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3 Results and discussion 

3.1 Bulk YDC conductivity 

Figure p] presents the activation energies, E a (eV), for vacancy migration in 
bulk YDC obtained from DIPPIM simulations using both the current pa- 
rameters (Table [I]) and those from our previous work [55], as well as, from 
experimental data available in the literature. The E a values predicted by the 
h-DIPPIM potential were calculated from Arrhenius plots of the ionic conduc- 
tivity in the temperature range 1273 K - 1673 K (NVT ensemble) for each of 
the three concentrations considered in this study. The previous LDA-DIPPIM 
(1-DIPPIM) work used 4x4x4 cells under NPT conditions [55J and a wider 
temperature range (1073 K - 2073 K). As expected, both sets of computational 
data are very similar despite being parameterized using different DFT func- 
tionals (h-DIPPIM shown as green circles and 1-DIPPIM as black triangles). 
Moreover, the values from the simulations are in good agreement with the 
range of available experimental results (blue diamonds |75j and red squares 
|76|). Table M contains the DIPPIM calculated bulk Ui (Scm -1 ) as a function 
of dopant concentration at 1473 K. Both the h-DIPPIM and the 1-DIPPIM 
potentials predict x = 0.12 to be the ionic conductivity maximum in YDC. 
The same level of agreement was observed between the potentials for all tem- 
peratues studied. This indicates that the h-DIPPIM potential maintains the 
previous potential's ability to account for the properties of YDC, but with 
the added advantages discussed above (Section J2J). As a result, all subsequent 
parts of this study were performed using the h-DIPPIM potential. 



x in Cei_ :c Y a ;02_^/2 


1-DIPPIM 


Oi (Scm- 1 ) 
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0.136 
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Table 3 Bulk ionic conductivities, cr, (Scm" 1 ), at 1473 K from LDA (1-DIPPIM) [581116] 
and HSE (h-DIPPIM) DIPPIM potentials (present work). 



Application of isotropic strain to bulk YDC caused an increase in the ma- 
terial's ionic conductivity, Oi. Figure [3| presents Oi of bulk Cei_ :E Y :E 02_ a ;/2 
when e = 0.021 relative to e = 0.000 (ct 21 /ct°) as a function of temperature. 
The conductivity maximum remains at x — 0.12 (red squares) for all the tem- 
peratures considered, however the highest increase in ai with strain (3.50 x) 
is observed for x = 0.18 (blue triangles) at 1273 K. These results also show 
that the impact of strain on Oi increases as the temperature decreases from 
1673 K to 1273 K. This effect can be explained in terms of the decrease in 
the activation energy of vacancy migration, E a (eV) (Table El) due to the ac- 
tivated nature of the conduction mechanism. The magnitude of the changes 
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Fig. 2 Bulk activation energies, E a (eV), calculated using the h-DIPPIM (green circles) 
and 1-DIPPIM potentials (black triangles. Experimental bulk activation energies from Ou 
et al. I75II77I (blue diamonds) and Wang et al. |76| (red squares) 



observed in E a for YDC is similar to those reported by De Souza et al. [57], 
who found a decrease of up to 40% for CeC>2 at e ~ 0.02. Lower temperatures 
were not included in this study as they are more difficult to simulate because 
they require longer MD trajectories (in the order of several ns of simulation 
time). Unfortunately, isotropic strain is difficult to realize experimentally and 
is presented here due to its theoretical interest. As a consequence, the changes 
in conductivity experienced by YDC under these conditions will be used only 
as a benchmark when compared to those of anisotropically strained slabs in 
the next section. 
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in Ce 1 - x Y x 2 _ x/2 
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0.007 


E a (eV) for e 
0.014 


0.021 


0.08 
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0.487 [-14.6] 


0.397 [-30.4] 


0.338 [-40.7] 


0.12 
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0.554 [-14.8] 


0.392 [-39.6] 


0.18 
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0.567 [-15.9] 


0.493 [-26.8] 



Table 4 Activation energies, E a (eV), as a function of strain, e, for bulk YDC. The values 
in square brackets correspond to the percentage change with respect to the activation energy 
without strain for a given concentration. 
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Fig. 3 Ionic conductivities of bulk Cei_ x Y :E 02_ a ;/2 when e = 0.021 relative to e = 0.000 
(<7 2,1 /cr°) as a function of temperature. The different concentrations are indicated as x = 
0.08 (black dots), 0.12 (red squares) and 0.18 (blue triangles). 



3.2 YDC slabs: surface effects 



The YDC slabs discussed in this section serve as models for the study of ionic 
conductivity of epitaxially grown thin films, in particular, coherently grown 
films where the lattice mismatch does not exceed e = 0.03. Figure [4] presents 
the 2D ionic conductivities, Oi (Scm -1 ), at 1273 K for unstrained YDC as a 
function of dopant concentration for bulk (blue triangles) and slabs (black di- 
amonds). These conductivities were calculated from the diffusion coefficients 
along the x and y directions only (D xy - 2D lateral diffusion) . The MD simula- 
tions reveal that, for zero strain, the conductivities of the slabs are 50% higher 
than those of bulk YDC of the same concentration. This pattern remains un- 
affected when all three components of the diffusion coefficient (D x , D y , D 2 ) 
are included in the calculation of the conductivity of bulk YDC and when this 
value is obtained from the average of the conductivities along the xy, yz, xz 
planes (3D conductivity shown as red dots). To the authors' knowledge, this 
is the first time that such a behaviour is observed in YDC. Apart from the 
increase in the conductivity, the behaviour of the two systems is similar, with 
the highest conductivity achieved for a dopant concentration of x = 0.12. 

The MD technique makes it possible to perform a detailed examination of 
the conductivity within each atomic layer along the 2-axis of the slabs (de- 
limited by horizontal blue dashed lines on the left in Figure [5]). The diffusion 
coefficients for each of these 13 layers of Ceo.8sYo.12O1.94 at 1273 K are pre- 
sented on the right hand side of Figure [5] as a function of depth (black dots 
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Fig. 4 2D (lateral) ionic conductivities, tn (Scm -1 ), for unstrained bulk (blue triangles) 
and slabs (black diamonds) as a function of dopant concentration at 1273 K. Bulk 3D ionic 
conductivities are shown as red dots. 



aligned to their corresponding layers). It is immediately clear that the layers 
with exposed surfaces show a significantly higher conductivity (~4x) than 
the bulk region of the slab, with the lowest conductivities seen in the subsur- 
face layers. The diffusion coefficient for bulk YDC of the same concentration 
as the slab is shown for comparison (vertical red dashed line). The increased 
conductivity in the surface regions is likely due to the more disordered atomic 
configurations observed in these layers. This higher degree of disorder is evi- 
dent when the positions of the oxide anions on the surface regions (solid blue 
circle in Figure pi) are compared to those in the low conductivity subsurface 
(dashed black circle in Figure [5]). 



3.3 YDC slabs: strain effects 



All combinations of strains, where e x ,e y — {0.007,0.014,0.021}, were ap- 
plied to the YDC slabs. However, it was found that, although the conduc- 
tivity tended to increase whenever there was strain along either direction, the 
change was generally larger when e x — e y (which is to be expected since YDC 
adopts a cubic fluorite structure). For this reason, we will only report strain 
where e x = e y for the slabs. In Figure |6| we report the conductivity vs strain 
for Ce 1 - x Y x O< 2 _ x j 2 slabs at 1273 K, where x = 0.08 (black dots), 0.12 (red 
squares) and 0.18 (blue triangles). The plots show that er^ increases with strain 
for all concentrations of YDC, and just as in the bulk system, x — 0.12 dis- 
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Fig. 5 3D diffusion coefficients for each atomic layer along the z-axis for Ceo.8sYo.12O1.g4 
at 1273 K (black dots). The vertical red dashed line indicates the bulk conductivity at 
the same dopant concentration and temperature. Attention is drawn to the high degree of 
disorder in anion sublattice on surface layers (solid blue circle) compared to the ordered 
subsurface layers (dashed black circle). 



plays the highest conductivity. The inset shows the change in conductivity of 
the same YDC systems when e x = e y = 0.021 with respect to the unstrained 
slabs (tr 21 /cr°) for all temperatures under study. As was seen in bulk YDC 
(Figure J3J), applying strain to YDC causes a larger increase in conductivity at 
lower temperatures. In the case of slabs, this change corresponds to a 1.44 x 
rise with respect to the unstrained slab at the same temperature for x = 0.18 
(c/. the increase for isotropically strained bulk was 3.5 x). This increase may 
appear modest when compared to the orders of magnitude seen in some of 
the YSZ literature [37], but it is in very good agreement with the latest re- 
ports on ceria and zirconia-based electrolytes 27,28,56,38 . Our findings are 
in also reasonable agreement with previous calculations 57.461147]. when the 
differences in strain values and temperatures are factored in. The activation 
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energies, E a (eV), of the YDC slabs (Table [5]) decrease as the strain increases, 
nevertheless, these changes are not as pronounced as was the case of bulk YDC 
(TablelS, especially when e = 0.021. Again, since the ionic conductivity is an 
activated process, the conductivity enhancement is expected to be higher at 
lower temperatures, as observed by Kushima and Yildiz for YSZ [H)ll47| . 
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Fig. 6 Ionic conductivities of Cei—^Y^C^- 
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slabs at 1273 K as a function of strain. The 



(blue triangles). Inset: 2D ionic conductivities of Cei—^Y^C^- 



0.08 (black dots), 0.12 (red squares) and 0.18 



x/2 



slabs when 



0.021 relative to e x 



0.000 (cr 2 - 1 /<r°) as a function of temperature. 



x in Cei_ a: Y ;c 02_ a; /2 


0.000 


0.007 


E a (eV) for e x = e y 
0.014 


0.021 


0.08 


0.511 


0.479 [-6.3] 


0.477 [-6.6] 


0.419 [-18.1] 


0.12 


0.608 


0.528 [-13.2] 


0.510 [-16.1] 


0.504 [-17.1] 


0.18 


0.691 


0.643 [-6.9] 


0.588 [-14.9] 


0.499 [-27.7] 



Table 5 Activation energies, E a (eV), as a function of strain, e, for YDC slabs. The values 
in square brackets correspond to the percentage change with respect to the activation energy 
without strain for a given concentration. 



The ionic conductivities of the YDC slabs manifest a high degree of anisotropy, 
which is evidenced when the diffusion coefficients , D (10~ 8 cm 2 /s), are de- 
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composed into the individual components along each Cartesian direction. The 
values reported in Table [6] correspond to the diffusion coefficients along the 
x (Dj.), y (D y ) and z (D z ) directions for biaxially strained Ceo.8sYo.12O1.94 
slabs and isotropically strained bulk (square brackets) of the same concentra- 
tion. The results show that D^ and D y increase more slowly in the slabs than 
they do in bulk when strain is applied. This anisotropy likely arises from the 
compression along the z direction of the slabs (Poisson effect). The D z values 
of our simulation slabs are constrained by the non-periodicity along this direc- 
tion and they will ultimately tend to zero as a result, provided the simulation 
time is long enough. However, in the timescales presented here the results in- 
dicate that the relaxation perpendicular to the surface plane does not lead to 
a diffusion enhancement. 



€x — £y 


D* (10-* cmVs) 


D y (10~ s crrrVs) 


Dz (10" s cmVs) 


0.000 


33.5 [23.0] 


33.7 [26.8] 


23.3 [24.4] 


0.007 


38.1 [42.9] 


37.5 [42.4] 


24.5 [39.9] 


0.014 


40.5 [59.2] 


44.1 [62.5] 


25.4 [59.1] 


0.021 


49.3 [89.8] 


47.1 [89.9] 


26.7 [87.0] 



Table 6 Slab diffusion coefficients along each Cartesian axis for Ceo.8gYo.12O1.94 at 1273 K 
under different strain, e, levels. Isotropically strained bulk diffusion coefficients are shown 
in square brackets for comparison. 



4 Conclusions 



This work attempts to rationalize the impact of strain in epitaxially strained 
RE-doped ceria. The simulations used accurate dipole-polarizable interatomic 
potentials derived directly from ab initio calculations. The simulation condi- 
tions were realistic from the point of view of the strains applied, dopant cation 
concentrations and the temperatures employed. Also, for the first time, the slab 
method was applied to this problem which allowed for a more realistic repre- 
sentation of biaxially strained systems. In addition, the results from the slab 
calculations illustrated the effects of surfaces on the ionic conductivity. The 
results obtained from our calculations indicate that there is a clear enhance- 
ment of the conductivity in the surface region of a thin film. This resulted in 
our films being ~ 50% more conducting than the corresponding bulk sample, 
for zero strain. This effect is related to the more disordered atomic configu- 
rations observed in the surface layers and the presence of under-coordinated 
atoms. When strain is applied, either isotropically (bulk) or biaxially (slab), 
a moderate enhancement in the ionic conductivity is observed, for the stud- 
ied temperature and strain ranges. Biaxially strained thin films are found to 
be less conductive than isotropically strained YDC, which has been ascribed 
to the highly anisotropic diffusion mechanism observed in the thin films. Our 
findings confirm those from recent experimental studies which employed films 
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of high crystallographic quality and found only a limited impact on the ionic 
conductivity from tensile strain [281138]. 
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